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INTRODUCTION 

A question which invariably arises when one considers the calculation of turbu- 
lent shear flows is, "How complex a model should be used to calculate such motions?" 
Available at the present time are models varying in complexity from very simple 
eddy -transport models to models in which all the equations for the nonzero second- 
order correlations are solved simultaneously with the equations for the mean vari- 
ables. For this reason, it might be instructive to present a discussion of the rela- 
tionship between these two models of turbulent shear flow. Two types of motion will 
be discussed: first, turbulent shear flow in a stratified medium and, second, the 
motion in a turbulent line vortex. These two cases are instructive because in the 
first example eddy -transport methods have proven reasonably effective, whereas in 
the second, they have led to erroneous conclusions. 

It is not generally appreciated that the simplest form of eddy -transport theory 
can be derived from second-order closure models of turbulent flow by a suitably lim- 
iting process. This paper will discuss this limiting process and the suitability of 
eddy -transport modeling for stratified media and line vortices. 

SYMBOLS 


a,b 


model parameters 


specific heat at constant pressure 


D 


operator 


g 


gravitational acceleration 


gpgk 


general acceleration vectors 
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i,j,k 

m,n 

N 

N Ri 

P 

P 

Q 

q 

r,<M 

r c 

T 

t 


indices 

proportionality constants 
stability parameter (eq. (90)) 
Richardson number 
parameter (eq. (49)) 
pressure 

turbulent energy, (UU + W + WW)^^ 

r • » 

scalar velocity, (ui)'uy 
cylindrical coordinates 
vortex core radius (fig. 7) 
temperature 
time 


U,V,W nondimensional second-order velocity correlations 

U,v,w mean velocity in r-, <p-, z-direction, respectively, for two-dimensional line 

vortex; mean velocity in x-, y-, z-direction, respectively, for atmospheric 
motion 


x,y,z Cartesian coordinates 

a length scale proportionality constant 

r vortex strength 

A difference 

5 breadth of layer under consideration 
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Kronecker delta 


6 ik 

e local deformation in vortex 

k Von Kar man's constant 

Ai,A 2,A3 length scales 
X dissipative scale 

fj. viscosity 

v kinematic viscosity 

p density 

Ty stress tensor 

T t turbulent shear stress 

Subscripts: 

char characteristic 

crit critical 

max maximum 

o undisturbed, adiabatic atmosphere 

Bars over a quantity indicate mean values. Primes indicate the instantaneous 
fluctuation of the quantity from its mean value. 

A SECOND-ORDER CLOSURE MODEL FOR TURBULENT SHEAR FLOW 

In reference 1 the author presented a discussion of the development of an invariant 
second-order closure model for turbulent shear flow in an incompressible medium. The 
basic equations of this model are 
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w + u '\r-p p ,i + L'i' " vuv u n,j 


— + = ” u j u k 5f - ^ Sjf + % A 2 q (fxj u i u k + w { u j u k 


+ u j u i) + fed A 3 q ) + l^l A 3 q 


9 u j u i \ . a 


aujiik 




f ij = ^( u U +u j,i) 


q = |(u l ) , uj 


In this model, the length scales Aj, A 2 , and A 3 are assumed to be proportional to one 
another. The dissipative scale X is given by 

A 2 

x 2 = ( 6 ) 

a + (bpqAj/p) 

From a rather lengthy parameter search to determine the values of the quantities a, b, 
Ai, A 2 , and A 3 that are to be used in the calculations, it was found that good results 
were obtained by using 

a = 2.5 " 

b = 0.125 


E 


It was further found that the remaining free parameter Aj was approximately equal to 
0.6 times the longitudinal integral scale of the motions studied in the parameter search. 

In reference 2, this model was extended to the case of turbulent motion and trans- 
port in the earth’s boundary layer. The case considered is that when the scales of the 
mean distributions of velocity and temperature are not greatly different and the Prandtl 
and Schmidt numbers are 1. The resulting equations, written in Cartesian tensor nota- 
tion, are 




In these equations, 


D _ 

DF" 


aF + u j a£ 


u . 9 

j 


(14) 


T 0 and v Q are the local temperature and kinematic viscosity, respectively, in an undis- 
turbed adiabatic atmosphere; T is the departure of the mean temperature from the 
adiabatic temperature T 0 ; and T' is the instantaneous fluctuation of the temperature 
about its mean. 


For the two cases of motion considered herein, equations (1) to (14) reduce to the 
following: 

For the case of a two-dimensional line vortex in cylindrical coordinates r,$,z 
with velocities u,v,w, equations (2) and (3) result in 


|f = -ffiV -JuV + i/te + Ig-^ 

at 3r r l ar 2 r ar r I 
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(18) 


aw'w* _ d_ [ a ~ dw'w' 
at " ar \ A 2 q ar 


. A 2 q aw’w’ 
+ ~d¥ 




a^w’w’ 1 aw'w' 2wV\ 
dr 2 rar ~~J~J 


auV 

at 



- 

ar 


+ 2 



auV" ) 2 ^Agqu’w’) 4A 2 q ^ 
3r / r ar r ar 




2 9(A 3 qu’v’) 
+ r ar 


A 3 q 8uV 
r ar 


4Agqu’v’ 


- uV + v 
A 1 


a 2 uW 

ar 2 


1 auV 
+ r ar 


4u'v' 

r 2 



(19) 


For the atmospheric motion considered herein, in which only a mean velocity u in 
the horizontal direction x exists and in which the mean lateral and vertical velocities 
v and w in the directions y and z, respectively, are zero, the appropriate equations 
derived from equations (8) to (13) are 


au _ ap .. a 2 u auV 

af ax o dz 2 dz 


( 20 ) 


aT _ a 2 T aT’w' 
at “ o 8z 2 az 


( 21 ) 
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~ 9u . g rTTfT , « 9 /» _ 9u’wM + la_/ A au'w' 
w az + T^ u T +Z 9 z l A 2 q 9z I P 0 9z l p o A 3 q 


2_t„ * 


.iuV + . 


0 az 2 


° x 2 


(25) 


au'T’ _ : 
9t 


uV f- w,T 'I + fe A 2 q az 


au’T’ 




Al 


az' 


(26) 


aw*T’ 

at 


- wV If + ’ ft< T ’> 2 + 2 ( A 2i It 11 ) ♦ h M p ° a 3“ az 


aw’T' 


q ir’T’ i i- ^ T *>,. W T 

T“ W 1 + t'o s iSl'o — 9— 

A 1 0 az 2 x 2 


(27) 


at az az 




9 2 (T*) 2 o (T’) 2 

+ — s — - — 2v„ - — — 


0 az 2 


0 X 2 


(28) 


For most people who practice the art of predicting turbulent flows, the equations 
for the line vortex (eqs. (15) to (19)) and for the flow in the atmospheric boundary layer 
(eqs. (20) to (28)) are of a familiar form. If distributions on the dependent variables are 
given at time t = 0, the development of the motion at subsequent times can be computed 
by simultaneous solution of the appropriate coupled sets of partial differential equations. 
This is the general approach of second-order modeling. 

The older and still widely used method of treating these problems is to consider 
only the equations for the mean variables and to assume that the second-order correla- 
tions which appear in them might be represented by empirically determined eddy -transport 
models patterned after the transport of the appropriate quantity by molecular means. The 
section which follows will examine what information can be gleaned from the equations for 
the second -order correlations about the nature of such eddy -transport models. 

SUPEREQUILIBRIUM MODELS 

To determine how to obtain information about the nature of eddy -transport models 
from the model or rate equations for the appropriate second-order correlations, one must 
consider what is implied when it is assumed that a turbulent flow can exhibit an eddy vis- 
cosity or an eddy diffusivity. 
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First, it is apparent that if the turbulent transport of a quantity depends only on the 
local gradient of that quantity and a scale length associated with the mean flows at the 
location under consideration, the turbulent transport cannot have a "memory" of its past 
history along the streamline. This is tantamount to the assumption that at each point in 
the flow the turbulent -transport correlations can track their local equilibrium values. 
These local equilibrium values can be obtained from the rate equations for the correla- 
tions by setting the left-hand sides of the equations, as they are given in the preceding 
section, equal to zero. Thus it is assumed that the rate of change of a transport correla- 
tion as it follows the mean motion is small compared with the production, dissipation, and 
diffusion terms which occur at the point in question. 

Second, the notion of an eddy -transport coefficient is one which does not allow the 
behavior of the turbulent transport at one point in the flow to affect directly the turbulent 
transport at another point. This notion is equivalent to the neglect of the diffusion terms 
in the equations for the second-order correlations, for it is these terms which link the 
generation of transport correlations at one point in the flow to the transport correlations 
at another point. 

Finally, the use of an eddy -transport model is a practice generally restricted to 
flows with high Reynolds numbers. Therefore, the high Reynolds number limit of the 
equations for the second-order correlations can be taken if it is desired to derive a sim- 
ple form of eddy -transport model hom these equations. 

If the three rules set forth above are followed, it should be possible to derive from 
the equations for the second -order correlations a simple theory of eddy transport. As 
discussed above, this theory represents the equilibrium, nondiffusive, high Reynolds num- 
ber limit of a second -order closure model. For reasons of brevity, this limit has for 
some time been referred to by the author as the "superequilibrium" limit. 

By following the three rules set forth, the following equations are found to be the 
superequilibrium equations for a line vortex: 


0 = ^uV - (1 +2b)SLuV +|_ 


0 = 


0 = 


0 = 




-(1 + 2b) ^ w'w’ + 7 ^ 


- fjp + jju’u' + ^ v'v' - (1 + 2b) % uV 


(29) 

(30) 

(31) 

(32) 
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For atmospheric motion, one finds 


3 

-2u'w' - 1(1 + 2b)uV + 1- 

9z A v ' A 

(33) 

-1(1 +2b)?V 

(34) 

ijp-w’T’ -1(1 +2b)wV 

(35) 

-wV + J- u’T’ - 1(1 + 2b)u’w’ 

o Z 1 q A 

(36) 

•uV |1 - w'T' |M _ 1(1 + 2b)u'T' 
9z 3z A 

(37) 

-w V |1 + ^(T’) 2 - 1(1 + 2b)w’T' 

(38) 

-2w’T’ |1 - 2b 1(T') 2 
3z A 

(39) 


In writing these equations, Aj was taken equal to A. 

EDDY TRANSPORT IN THE ATMOSPHERE 


It is instructive to carry out the solution of equations (33) to (39). These equations 
are algebraic for all the nonzero correlations. The solution of equations (33) to (39) can 
be obtained if the following definitions are introduced: Let (for 9u/ 3z > 0) 


o/ -\ 2 
u’u' = UUAj 2 ^) 

u’T' = UTAj 2 f |1 

Tv = wai 2 (||) 2 


2 

w V = WWA X 2 ) 

(T’) 2 = TTAj 2 

_ 2 

uV = UWA! 2 (||) 

q 2 = QQA! 2 (f) 2 


(40) 


and note that 


QQ = Q 2 = UU + VV + WW 


(41) 
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Substitution of the definitions given in equations (40) and (41) into equations (33) to (39) 
results in 


Q(1 + 2b)UU = &- - 2UW 

(42) 

C)3 

Q(1 + 2b)VV = 

u 

(43) 

Q(1 + 2b)UW = -WW + N Ri UT 

(44) 

Q(1 + 2b)UT = -UW - WT 

(45) 

Q(1 + 2b)WT = -WW + N Ri TT 

(46) 

Q(2b)TT = -WT 

(47) 


In these equations, N Ri is the Richardson number given by 


N Ri = 


_g_ 8T 
T 0 3z 



(48) 


It is immediately obvious from equations (42) to (47) that all the nondimensional second- 
order correlations are a function only of the Richardson number and the parameter b 
from the second-order closure model. It will be remembered that the value of b deter- 
mined in the parameter search reported previously is 0. 125. 

It is convenient to express the solution of equations (42) to (47) in terms of the 
parameter 


P = 


1 - (4 + 15b)N m + 


1 + 2(2 - 9b)N Ri + (4 + 9b)' 


Kif] 


1/2 


(49) 


In terms of this parameter, the various correlations may be written 



b(l + 2b; 


(50) 
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All the nondimensional second -order correlations as functions of the Richardson 
number are plotted in figures 1 to 5. From these figures, the profound difference between 
turbulence and turbulent transport in stable and unstable atmospheres is obvious. Note 
particularly that the nondimensional vertical transport of matter and heat falls off far 
more rapidly than do the nondimensional turbulent energy components when a stable 
atmospheric situation is approached. In fact, above a Richardson number of 1, vertical 
turbulent transport has almost ceased to exist although there is still some atmospheric 
turbulence. 
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It should be noted that the superequilibrium results just obtained specify the non- 
dimensional values of second-order correlations. For example, if the value of b = 0. 125 
is assumed to be correct, a Richardson number of 0.10 would give 


UW = -0.2712 

(60) 

WT = -0.2631 

(61) 


The transports of momentum and heat would then be given by (for du/dz > 0) 


-p o uV 


= 0.2712p o A 2 



(62) 



0.2631p Q A 


2 9u 9T 

dZ 3Z 


(63) 


It is clear from these expressions that the actual transport is not defined until the length 
scale A is known. This is a difficulty with atmospheric flows, for unless A is deter- 
mined at a given altitude and the local Richardson number specified there, the transports 
are not known. In general, A will depend at a given altitude on the Richardson number 
but can assume a range of values depending on the past history of the motion. Although 
this range of values is limited so that the order of magnitude of the transport might be 
determined, there will always be a variation in transport proportional to the square of the 
variation in A at any fixed Richardson number. 

For classical laboratory flows, this problem does not exist. In this case, it is gen- 
erally found that A is proportional to the characteristic breadth of the layer under con- 
sideration while the gradients are proportional to a characteristic velocity, temperature, 
or concentration difference divided by this characteristic breadth. Thus, for the classi- 
cal shear flows, 

Al 2 (f f - 6 char = Const (A“char) 2 < 64 > 


and, likewise, 


A l 2 H fz" “ Const Au char AT char 

For each type of flow, these constants are well defined. This type of simplicity is, alas, 
not true of the atmosphere. 
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It is instructive to compare the results of superequilibrium theory with certain 
well-known results from classical turbulent -transport theory for the case when no grav- 
itational effects are involved. To do this, the Richardson number is placed equal to zero 
in the expressions given in equations (49) to (57), and for b = 0.125, the following expres 
sions are obtained: 


P = 0.3333 


Q 2 = = 1.7066 

3b(l + 2b) 2 

UU = — 1 + 6b = 0.7964 
9b(l + 2b) 2 

VV = WW = 1 = 0.4551 

9b(l + 2b) 3 

UW = TW i ffik = -0.2786 

9(1 + 2b) 3 

UT = 2 =•= 0.3413 

3(1 + 2b) 3 

TT = 1 = 1.7066 

3b(l + 2b) 2 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 


Some interesting results are n oted f rom the above comparison. First, superequi- 
librium theory indicates that v'v' = w'w' and, further, that 



v'v' 



w'w' 


UU 

VV 


UU 

WW 


= 1 + 6b = 1.75 


(73) 


Second, the value of -u'w'/q 2 , which Bradshaw, Ferriss, and Atwell (ref. 3) assume to 
be a constant equal to 0.15, is defined by superequilibrium theory to be 


This is a rather surprisingly accurate result in view of the fact that the value of b was 

determined from very different considerations in the development of the second-order 
closure model. 
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The value of Von Kar man's constant k in his expression for the turbulent shear 
near a surface may also be derived from superequilibrium theory: 


T t = 


■pu'w' = PK 


2 Z 2 



(75) 


From the present results, 


T t = -pUWA^ff) 2 


\/3?b 

9(1 + 2b) 3 



(76) 


In reference 1 it was found that near a surface, A is of the form A = az, where a in 
the parameter search is found to be 0.7. Letting 


A = 0.7z 


(77) 


equation (76) gives 


T t = 


0.491/3/5 2 

pz, 

9(1 + 2b) 3 



Comparison of equations (75) and (78) reveals that 


k 2 , °- 49 \/3^ = Q is? 
9(1 + 2b) 3 


(78) 


(79) 


or 


The value of Von Karman’s constant is actually 0.4. Again, the agreement between results 
obtained by taking the equilibrium, nondiffusive limit of the present second-order closure 
model of turbulent shear flow and the classical mixing -length theory is rather remarkable. 

EDDY TRANSPORT IN A VORTEX? 

If a scheme such as that pursued in the previous section for the superequilibrium 
equation for a line vortex is followed, the following definitions are introduced into equa- 
tions (29) to (32): 
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(81) 


- _x2 


u V = UUA 2 (g - l) 


v V = WA 2 (g - ¥ 
\dr r 


w V = WWA 


2/jjv v'* 
3r 


uV = UVA 2 


9v v 
9r " r 


-?) 

/8v v\ 

Var ■ rj 

q = QA " 

From this substitution, the following equations are obtained: 


9v v 
8r " r 


(82) 

(83) 

(84) 

(85) 


n 3 

Q(1 + 2b)UU = ^y + 4WN 

(86) 

q3 

Q(1 + 2b)W = - 2UV - 4UVN 

(87) 

n3 

Q(1 + 2b)WW = 

0 

(88) 

Q(1 + 2b)UV = -UU + 2(VV - UU)N 

(89) 

In these equations N is a stability number defined as 


N = 7/r 

dV V 

dr ~ r 

(90) 

The solution of equations (86) to (89) in terms of the parameters 

b and N is 

UU = Q 2 A - 4bN\ 

1 + 2b ^3 J 

(91) 

VV= lt2b(5 + 2b+4bN ) 

(92) 

ww = — 

3(1 + 2b) 

(93) 
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(94) 


UV = -bQ 3 

with 

Q 2 = UU + W + WW J\ - 8bN - 16bN 2 ^ (95) 

b(l + 2b) 2 v / 

It is clear, since Q 2 is positive definite, that under the assumptions made here 
turbulence is impossible if 

N< -i + (96) 

or if 

n> ?(iF^- 1 ) (97> 

For the value of b used herein (0.125), these limits become N< -0.729 and N> 0.229. 

Figure 6 shows the behavior of the quantities UU, VV, WW, UV, and Q 2 with 
variations of the stability parameter N for b = 0.125. The results are plotted in terms 
of the ratios of the quantities to their values for N = 0, nam ely, (UU) 0 , (W) 0 , and so 

forth. Thus, figure 6 shows the ratios of uV, v'v', w V, u'v', and q 2 in a vortex 
to these quantities in a parallel shearing motion having the same mean deformation rate 

and scale. 

It may be seen from figure 6 that the turbulent energy and shear have the same 
value for N = -1/2 as they do for N = 0. Between N=-l/2 and N = 0, the turbulent 
energy and shear are larger than they are in a parallel shearing motion. For N < -0.729 
and N > 0.229, as mentioned previously, no locally sustained turbulent flow is possible. 
Thus, for -0.729 < N < -0.5, locally self -sustained turbulence is possible, although the 
turbulence is damped by centrifugal effects. For 0 < N < 0.229, turbulence is also pos- 
sible, but here again it is damped by the action of centrifugal forces. 

What sort of flows does each of these regions represent? First, note that when 
9v/9r = 0, N = -1. Thus at the core radius of a vortex (defined here as the radius where 
9v/9r = 0), a turbulent vortex is stable. Near the center of a free vortex, the tangential 
velocity v is of the form v = mr - 2nr 2 so that as r - 0, N - -°°. Also, for a free 
vortex, v - r/27rr as r - « and one finds then that as r - N — -1/2. Thus for 
the classical vortex distribution, 
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(98) 


The flow in the outer regions of the vortex exhibits an eddy diffusivity similar to a par- 
allel flow. As the core of the vortex is approached, the flow becomes more and more 
stable. It becomes completely stable somewhat outside the core of the vortex. Indeed, 
the flow is stable at the point of maximum deformation |p - p This behavior of the 
stability parameter N for the classical vortex is shown in figure 7. 

The region of increased turbulence and shear between N=-l/2 and N = 0 can 

be understood if it is noted that the stability parameter may be written 
xt - v/r 

_1_ il _ 2v (") 

2tt 3r r 

Thus the region -1/2 < N < 0 represents flows for which dr/dr is negative. These 
are, of course, flows which exhibit the well-known Taylor instability (ref. 4). 

The region 0 < N < 0.229 is representative of flows occurring between two cylin- 
ders rotating in the same direction, so that r at the outer cylinder is larger thqn r 

at the inner cylinder when the centrifugal forces due to the general level rotation cannot 
completely stabilize the flow. 

For a free vortex, it may be surmised from this analysis that the core regions of 
vortices are locally stable. Regions outside the core are unstable and can generate tur- 
bulence. If the core regions of vortices are to exhibit a turbulent shear, this must be 
caused by turbulence which has diffused into the core region from outer regions which 
are unstable or by turbulence which has been generated by a shear in the axial direction 
that is not considered in this analysis. This fact, namely, that the turbulent shear -puV 
in a vortex is not directly related to the local deformation g - |, would lead one to 
believe that it would be impossible to establish any general rules for determining an eddy 
viscosity for a vortex. To calculate such flows reliably, it will probably be necessary 
to use the full power of second-order closure methods. 

It might be noted, in this connection, that if one were to use an energy method on 
such flows, much of the physics of the problem would be lost. This may be seen by con- 
sidering the sum of equations (86) to (89) with UU = VV as the governing equation of the 
flow. In this case, the parameter N disappears from the equations and the essential 
physics of the problem have been lost. 
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CONCLUDING REMARKS 


This short paper has tried to exhibit the relationship between a full second -order 
closure model for turbulent flow and the older eddy -viscosity models of such flows. It 
has been shown that classical eddy -transport theory can be obtained from a consideration 
of the equilibrium, nondiffusive, high Reynolds number limit of the equations of a second- 
order closure model. 

The nature of such limits has been discussed for the flow in sheared stratified 
media and for a line vortex. The limitations of an eddy -transport model of turbulence 
in a line vortex have been discussed through the use of the equations derived by the lim- 
iting process described herein. 
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N = -0.59l — V-- 

Ncrit <N<-0.5\ 
locally self- sustained 
but damped turbulence 


N < N cr j t no locally 
self- sustained 
turbulence possible 


€/€ max = (d v/dr-v/r)/(d v/dr - v/r) 

Figure 7. - Behavior of the stability parameter N for a classical 

line vortex (see eq. (98)). 





DISCUSSION 


G. L. Mellor : We just finished an experiment at Princeton a year or so ago which is a 
boundary layer on a flat wall followed by a curved wall and we did a complete set of tur- 
bulent measurements. It was very dramatic indeed; we could see the turbulence nearly 
shut off on the convex side of the wall. And then we did try to develop an expression for 
the eddy viscosity by balancing production with dissipation using these equations. You 
can, as you did, find a correction for the eddy viscosity in terms of the proper curvature 
parameter and it works very well. So I'll just put in two cents here that it seems to work 
well and it seems to compare with the rather dramatic measurements that we've devel- 
oped in the last couple of years. 

C. duP, Donaldson : One of the points to be made here is that if you are going to do some 
trick flow like this, it is best not to shortcut and just use, say, the energy and stress 
equations with some trick for guessing what is missing because you may overlook the 
physics of the problem. If you use all the equations, you get the physics right. 

S. I. Pai : I notice that you assume the mixing lengths for velocity and for temperature 
are the same. My questions are (1) Do you make this assumption to simplify the analysis 
and to obtain some essential features of this problem? and (2) If these mixing lengths are 
different, would you expect that your results would be modified considerably? 

C. duP. Donaldson ; Yes, the same lengths were chosen to simplify the analysis. In this 
case, which is the superequilibrium limit, you will begin to see the nature of the problem. 
No matter what you choose for your eddy -viscosity model, if you want to see whether you 
should really use such a model or not, you can make this kind of limiting argument. It is 
true that if you have vastly different scales of the temperature and velocity fields, you 
are going to have to do something different. As an example, take a large turbulent pipe 
flow in which you know all the pertinent mean quantities, and with the turbulence model 
you have, compute all the turbulent characteristics in this flow. Then, if you assume a 
tiny pencil of heated air is placed in the center of the tube so as to fo rm a very thin hot 
jet in that r egio n and you use the same scale in the u-T' and (T’) 2 equations as you 
use in the uju^ equation, you will find a remarkable result. The general spread of the 
hot material is about as it should be, but a hot spot stays ne ar the center of the tube. 

This is a result of using the wrong scale in the u'T’ and (T’) 2 equations. When the 
scales of the mean temperature and mean velocity fields are so disparate, one may not 
use the same scales in the equations for the various second-order correlations. 

S. C. Lee : I have two questions. First, I haven't seen any of your models compared with 
the suggested cases. My question is, have you compared any, and if you did, would we be 
able to see any of your comparisons? 
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C. duP. Donaldson : The only comparisons made to date were for the cases which were 
used to construct the model. In these cases, model parameter searches were made to 
try to obtain a best fit to the existing data for a free jet, a two-dimensional shear layer, 
and a flat -plate boundary layer. The model was determined to be that one that gave the 
best fit to all three cases. 

I think I should point out here something that I've said ever since I started making such 
calculations; namely, it's pretty hard to calculate better at first blush something that 
somebody has been calculating empirically for the last 20 years. That really wasn't my 
reason for getting into these second -order closure methods - my real reason was to get 
at the calculation of some problem which just can't be done by conventional methods and 
for which you just don't known what the answer is at all, such as the vortex and the behav- 
ior of turbulence in stratified media. 

S. C. Lee: My second question is related to your particular model with which you are 
interested in atmospheric conditions for stable and unstable atmospheric conditions. 

Have you calculated any of those? 

C. duP. Donaldson: Yes. 

S. C. Lee : Would those be in the paper? 

C. duP. Donaldson; Not in the paper to be published in these proceedings. Some results 
have been published elsewhere. I have just finished writing a paper which is to be mailed 
out soon to many of the people at this meeting. I think you are on the list. 

S. Corrsin : You have identified one necessary condition for the use of an eddy -viscosity 
model and you of course know that also there is another kind of necessary condition; that 
is, the characteristics length of the mechanism transporting the property you are inter- 
ested in must be very small compared with any distance over which the mean property 
changes appreciably and this is violated by almost all flows that we ever talk about. So 
there is quite a different class of reasons why eddy -viscosity models may be wrong in 
principle but work in practice. 

C. duP. Donaldson : I understand - valid comment. 

G. L. Mellor : I ask Stan, who states this necessary condition - I know it has something 
to do with kinetic theory but who states it for turbulent? 

S. Corrsin: Well as far as ordinary transport models go there's a book on Transport by 
a fellow named Bosworth* that was published in the 1940' s which is the only textbook I 
have seen that actually mentions it. But as far as turbulence phenomena go, in general the 
gradient transport term is the first term in an infinite -series approximation and this is 

* Bosworth, Richard Charles Leslie: Heat Transfer Phenomena - The Flow of Heat in 
Physical Systems. John Wiley & Sons, Inc., [l952]. 
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why, for instance, the Prandtl modified-mixing-length theory with the second-order term 
could be better. In fact, if your computers get big enough, it may be better to just put on 
higher and higher derivatives in an attempt to make something which is physically mean- 
ingful as well as computable. 

C. duP. Donaldson : Yes, that’s certainly true for those cases that you can do like that. 
When, indeed, the flow is completely stable in the superequilibrium sense in the region 
where the deformations are largest, you can't do that. 
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